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Deep optical surveys [j^, |2[] probe the rest—frame ultraviolet luminosities of 
high—redshift galaxies, which can be converted into star formation rates under 
plausible assumptions on young stellar populations. The current analysis of 
these data suggests that the global star formation rate of the universe peaked at 
a redshift of 1 and declined since then P, ^ . This has led to claims that the bulk 
of star formation in the universe has been seen. However, the conversion of UV 
luminosities into star formation rates must take into account a correction for 
the luminosity fraction absorbed by the dust which is generically associated to 
young stars. Since this correction is rather uncertain for high—redshift galaxies, 
the star formation rates currently deduced from optical surveys alone might be 
substantially underestimated. To circumvent this problem, the simplest is to ob¬ 
serve the dust thermal emission at infrared (and submillimetre) wavelengths and 
compute the overall luminosity bugdet of galaxies. For high—redshift galaxies. 
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the only direct observational constraint is set by the recent detection of the Cos¬ 
mic Infrared Backgronnd (CIRB) built up from the accumulated IR light of faint 
galaxies along the line of sight [§. Here we propose a more accurate determi¬ 
nation of this long—sought background which solves the main possible weakness 
of the earlier determination. Then we estimate the population of high—redshift, 
dust—enshrouded starburst galaxies needed to produce this background. We ar¬ 
gue that most of the star formation at high redshift may be hidden by dust, and 
we define the necessary characteristics of a feasible survey at a wavelength of 
175 /xm, which could detect this population. 

While only one third of the bolometric luminosity of local galaxies is radiated in the IR 
1^, there is a growing evidence that this fraction is actually increasing with redshift. The 
deepest counts available from the Infrared Astronomical Satellite (IRAS) at 60 /im [0], which 
correspond to an average redshift z of only 0.2 [^, already suggest some evolution of the IR 
emission in the universe. A recent deep survey with the Infrared Space Observatory (ISO) 
at 15 pm has discovered a few objects at z ~ 0.5 to 1 with star formation rates much higher 
than deduced from the optical P]. However, the strongest constraint on the high-redshift IR 
emission is given by the CIRB found in data acquired by the FIRAS instrument on-board 
the COBE satellite in the 200 pm - 2 mm wavelength range [Q . Several steps are necessary 
in order to remove the foregound Galactic components and extract the isotropic residual 
identified as the CIRB. While, at the wavelengths probed by FIRAS, the interplanetary 
emission is small and easily removed 0, the emission from interstellar dust mixed with 
the different gas phases of the interstellar medium is the dominant component. The spectra 
which correlate with the 21 cm interstellar emission of neutral hydrogen (HI) []T| along 
lines-of-sight with HI column densities N^i < 4.5 x 10^° atoms cm“^ are described by a 
modified black-body v‘^By{T) with T = 17.5 K |]^]. Part of the long-wavelength excess 
over this simple model increases with HI column density and can be linked to dust emission 
associated with ionized and molecular hydrogen. The other part is an isotropic residual which 
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has been interpreted as the CIRB. The determination of the isotropy required to use a large 
fraction of the sky. This thus involved correcting for substantial foreground components. 
Any inaccuracy in this correction might have contributed a spurious signal. In order to 
address this problem, the original method of Puget et al. 0 was applied again, but only in 
the cleanest regions with very low HI column densities {Nhi < 1 x atoms cm“^ instead 
of Nhi < 4.5 x 10^° atoms cm“^). In that case, the residual component totally dominates 
the emission (inset in hg.l). This demonstrates that it cannot be due to artifacts in the 
removal of interstellar emission. Fig.l also shows that the CIRB intensity per frequency 
decade = 7 10“® W m~^ sr“^ near 300 /xm is a factor of 5 higher than the no-evolution 
prediction obtained by a simple extrapolation of the IR luminosities of local galaxies. Its 
level is comparable to that of the “Cosmic Optical Background” estimated by summing up 
faint galaxy counts down to the deepest limit so far available which is given by the Hubble 
Deep Field [§. 

In order to break the background into the contributing sources, we need to model the 
IR/submm emission of galaxies. Galaxy formation and evolution can be briefly sketched 
as follows: Small fluctuations in the high-z universe grow by gravitational instability until 
they form dense clumps where the baryonic gas is collisionally heated. Where the density 
and temperature are appropriate, gas cools by emitting radiation, and the baryons pile up 
in cold cores whose hnal radii are set by angular-momentum conservation. Simultaneously, 
larger clumps form and encompass/accrete the previous generation of small clumps. Stars 
form in the cores and enrich the primordial gas via supernova ejecta. Part of the starlight is 
absorbed by various dust components which re-radiate at longer wavelengths according to 
characteristic IR spectra. The so-called “semi-analytic modelling” of these series of physical 
processes has been rather successful in reproducing the overall properties of galaxies in the 

We have elaborated an extension of this method to the IR/submm 


optical range 14 


range. Details of the modelling, and complementary predictions will be given elsewhere 
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Specifically, we assume a standard Cold Dark Matter cosmological model [Q with a 
Hubble constant Hq = 50 km s“^ Mpc“^, a density parameter Dq = 1, a cosmological con¬ 
stant A = 0, a baryon fraction = 0.05, and a normalisation = 0.67 for the power 
spectrum of linear fluctuations. The sensitivity of the semi-analytic modelling to cosmolog¬ 


ical parameters is known to be weak |^. In our reference scenario A, we consider a mix 
of two broad types of populations, one with low star formation rates (which reproduce the 
observational distribution of gas consumption time scales in disk galaxies the other 

proceeding in bursts with ten times larger star formation rates. We assume that bursts are 
triggered by interaction and merging of sub-galactic clumps [^], and increase with accord¬ 
ing to the fraction of pairs ^ This population of “mild starbursts” and “luminous 
UV/IR galaxies” (similar to “LIRGs” |^) dominate the optical background (fig.l). The 
IR luminosity-to-mass ratio is Ljr/M = ^Li^oiq/Mq for a typical starburst with = 0.33 
Gyr. The range of derived IR-to-blue luminosity ratios is characteristic of blue-band se¬ 
lected samples like the Ganada-France Redshift Survey (selected in the observer-frame 
Iab band, roughly corresponding to the B band at 2 ; ~ 1), or the high- 2 ; galaxies of the 
Hubble Deep Field. Fig.l also displays the predicted IR background, which is clearly barely 
compatible with the observed GIRB whose mean amplitude is twice higher. 

In order to assess how much star formation might be completely hidden by dust shrouds, 
we consider an additional population, similar to “ultra-luminous IR galaxies” (“ULIRGs”) 
25|. We maximize their IR luminosities by assuming that all the energy available from 


stellar nucleosynthesis is radiated by massive stars and heats up the dust. As a consequence 
of a stellar initial mass function with short-lived, massive stars, the post-starburst phase 
is “dark” and would be detectable only by its nucleosynthesis products. The luminosity- 
to-mass ratio now is Lj^/M = ISOL^oiq/Mq for a typical starburst with G = 0.33 Gyr. 
Our scenarios B and G respectively mimic continous bulge formation as the end-product of 
interaction and merging, and a strong episode of bulge formation at Zfor > 3.5. Both are 
consistent with the observed GIRB. The high-redshift, dust-enshrouded, star formation in 
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scenario C results in the high level of the predicted background at wavelengths A > 400 /iui. 

While none of the currently available optical data reflects the large differences between 
these scenarios, originating in the different fractions of heavily-extinguished objects, the 
predicted IR/submm counts are more interesting, as shown in hg.2. The comparison of 
IRAS data with the no-evolution curve at 60 /xm suggests some evolution. However, it 
appears that the 60 /rm band does not strongly discriminate between the various scenarios 
of evolution. In contrast, the upward deviation at 200 /im is due to the contribution of the 
redshifted 100 /xm maximum of the IR energy distribution. This redshifting of steep spectra 
counter-balances distance dimming and can make high-;^ objects easier to detect than low-;^ 
ones. Submm observations are thus quite sensitive to the high-;^ history. The model also 
predicts that, at 200 /rm, 10-100 mJy sources (contributing to 15 % of the background) are 
mostly located a.t z ^ 0.5 - 2.5, while at 60 /rm, and at the typical sensivity level of IRAS 
surveys, the sources are indeed located mostly at very low ; 2 . 

The detection of these sources would be a strong test for assessing the level of the 
“optically-dark” side of galaxy formation. The C160 hlter of the ISOPHOT instrument 
on-board ISO has an effective wavelength Ag// — 175 /rm for typical spectra of distant 
galaxies, and a 10 mJy rms noise fluctuation per 1.5 arcmin pixel is reachable after integra¬ 
tion times larger than ~ 256 s per pixel. Thus a deep survey with this instrument appears to 
be feasible and is indeed scheduled. However, one might be concerned that small-scale cirrus 
fluctuations could hide the sources and the fluctuations of the background they induce. A 
comparative analysis of the expected power spectra due to (1) cirrus fluctuations in regions 
of various HI column densities, (2) background fluctuations once sources above the confu¬ 
sion limit have been removed, (3) the detector noise, shows that, in clean regions of the sky 
{Nhi < 1 X 10^° atoms cm“^), a survey with 10 mJy rms sensitivity should not only detect 
most sources above the low cirrus fluctuations but could also see background fluctuations in 
excess of the detector noise fluctuations, on scales 3 to 10 arcmin (see hg.3). Since scenario 
C has 6.3 X 10^ sources/sr with fluxes > 30 mJy, a deep survey of a ~ 1000 arcmin^ field 
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might begin to “break” the CIRB into ~ 50 ± 7 discrete sources, an order of magnitude 
more than is expected without evolution. This number is sufficient to test the high level of 
evolution and even begin to disentangle between the various scenarios. If detected, the level 
of background fluctuations can also help to constrain the redshift distribution of the sources. 

The tentative discovery of the Cosmic Infrared Background is now supported by our 
new study in the cleanest regions of the sky, where the foregound Galactic components are 
essentially negligible. The conversion of the CIRB into its contributing sources by means 
of a semi-analytic model of galaxy formation leads to predictions of faint galaxy counts at 
IR and sub mm wavelengths. In contrast with the status of IRAS 60 pm counts, 175 pm 
counts with ISO should be able to detect the predicted strong evolution, and even disentangle 
between various scenarios consistent with the level of the CIRB. Moreover, small-scale cirrus 
fluctuations cannot hide the presence of the sources. Consequently, our knowledge of the 
optical/IR luminosity budget at z ^ 0.5 - 2.5 should improve rapidly. We are about to start 
unveiling the optically-dark side of galaxy formation. 
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Figure Captions 


Figure 1: Inset panel: High-latitude COBE/FIRAS spectrum in regions with very low 
HI column densities {N^i < 1 x 10^° atoms cm“^, solid line) and residual spectrum after 
subtraction of emissions correlated with neutral and ionized hydrogen (dotted line). This 
component is likely to be the Cosmic Infrared Background (CIRB). Main panel: The ±lcr 
error bars per point have been used to define an acceptable range for CIRB predictions 
(thick solid lines). Solid squares show the upper limits given by COBE/DIRBE residuals 
0. The solid hexagons show the Cosmic Optical Background obtained by summing up 
faint galaxy counts down to the Hubble Deep Field limit. Strictly speaking, this is only a 
lower limit of the actual optical background, but the shallowing of the U and R-band counts 
suggests near-convergence at least at those wavelengths p| . The short and long dashes show 
the background predicted without evolution, in a universe with Hq = 50 km s“^ Mpc“^ and 
Do = 1- The solid, dotted and dashed lines show the predictions for our scenarios A, B, and C 
respectively, in the standard CDM model (Rq = 50 km s“^ Mpc“^, Dq = 1, A = 0, D;, = 0.05, 
and (Ts = 0.67). Stars form in cold cores according to Salpeter’s stellar initial mass function 
(with slope s = 1.35, and lower and upper mass cut-offs 0.1 and 120 Mq). Star formation 
rates are proportional to the cold gas contents, with characteristic time scales derived from 
the core dynamical times as = jStdyn- Scenario A has a mix of two populations, a disk¬ 
like one with fl = 100, and starbursts with /? = 10. The mass fraction involved in bursts 
increases with the formation redshift Zfor according to /burst oc {1 + ZforY, as suggested by 
the increasing fraction of pairs seen at larger 0, |^, In scenario B and C, we add 
a population of “ultra-luminous IR galaxies” (ULIRGs) in which all the energy available 
from stellar nucleosynthesis (O.OOTxMc^) is radiated by massive stars (< x >= 0.40) in a 
heavily-extinguished medium. Scenario B has a 5 % constant mass fraction of ULIRGs at 
all Zfor- In scenario G, 90 % of all galaxies forming at Zfor > 3.5 are ULIRGs. 


Figure 2: Predictions for differential counts normalized to Euclidean counts at 60 pm 
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(upper panel) and 200 pm (lower panel). The short and long dashes show the predicted 
counts without evolution, in a universe with Hq = 50 km s“^ Mpc“^ and flo = 1- The solid, 
dotted and dashed lines show the predictions for our scenarios A, B, and C respectively 
in the SCDM model (see text). Data are plotted for IRAS counts at 60 pm. Open stars: 
Faint Source Survey Open squares: QMW survey ||2^. Solid squares: North Ecliptic 
Pole Region 0. The IRAS 60 pm counts suggest some evolution, but do not discriminate 
between the various scenarios. These scenarios also predict a 60 pm background fluctuation 
per beam in the Very Faint Source Survey (after removal of > dcxtoi = 120 mJy sources) 
at the level of 14.1 mJy (A), 16.0 mJy (B), and 14.3 mJy (C), while the measured 68 % 
quantile is 30.1± 1.2 mJy |^. With a 25 mJy rms instrumental noise and 6.5 mJy rms 
cirrus fluctuations |^, there is still space for a IbAtl’.l mJy fluctuation due to sources, in 
good agreement with our estimates. In contrast, scenarios A, B, C predict much stronger 
evolution at 200 pm. The differences in the predicted counts originate from the differences 
in the high-redshift IR emissions. 


Figure 3: Comparison of predicted power spectra for observations in the ISO C160 
hlter. The straight horizontal lines (same code as in hg.l) show the predicted background 
fluctuations at the confusion limit (i.e. its rms value equals one third of the limiting flux of all 
resolved and removed sources). The beam was modelled as a 1 arcmin FWHM gaussian. The 
cirrus fluctuations (thin lines) are described by a k~^ power law with levels (from left to 
right) set by column densities of neutral hydrogen (HI) respectively typical of the Lockman 
hole (minimum, 0.24 % of the sky), and of 4.6 %, & 21.8 % of the sky. The detector noise 
spectrum of fluctuations (thin dashed lines), assumed to be white, is shown for a sensitivity 
level of 10 and 20 mJy per 1.5 arcmin pixel. To allow for a direct comparison, this flat 
white noise has been projected on the sky, i.e. divided by the Fourier transform of the beam 
prohle. The resulting exponential rise of the on-sky noise thus bounds the range of available 
scales. The comparison shows that, in the cleanest regions of the sky, (e.g. Nhi = 1 x 10^*^ 
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atoms cm“^), unresolved background fluctuations for scenarios B and C actually dominate 
both cirrus and noise fluctuations at scales 3 < 0 < 10 arcmin, if the rms noise level is 10 
mJy per pixel. 
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